Mon. Not. R. Astron. Soc. 000, \T§E\ (2011) Printed 4 January 2012 (MN MgX style file v2.2) 



Broadband Observations of the Be/X-ray Binary Pulsar 
RX J0440. 9+4431: Discovery of a Cyclotron Absorption 
Line 



S. S. Tsygankov 1,2 ' 3 ' 4 *, R. A. Krivonos 3 ' 4 , A. A. Lutovinov 4 ' 2 

1 Finnish Centre for Astronomy with ESO (FINCA), University of Turku, Vaisdldntie 20, FI-21500 Piikkid, Finland 

2 Astronomy Division, Department of Physics, FI-90014 University of Oulu, Finland 

3 MPI for Astrophysics, Kari-Schwarzschild Str. 1, Garching, 85741, Germany 

4 Space Research Institute of the Russian Academy of Sciences, Profsoyuznaya Str. 84/32, Moscow 117997, Russia 



X 

6 



> 

q 
o 



Accepted 2011 December 30 



x 



ABSTRACT 

We report the results of an analysis of data obtained with the INTEGRAL. Swift 
and RXTE observatories during the 2010 April and September outbursts of the X-ray 
pulsar RX J0440. 9+4431. The temporal and spectral properties of the pulsar in a wide 
energy band (0.6 — 120 keV) were studied for the first time. We discovered a ~ 32 
keV cyclotron resonant scattering feature in the source spectrum, that allowed us to 
estimate the magnetic field strength of the neutron star as B ~ 3.2 x 10 12 G. The 
estimate of the magnetic field strength was confirmed by a comprehensive analysis of 
the noise power spectrum of the source. Based on the recurrence time between Type 
I outbursts the orbital period of the binary system can be estimated as ~ 155 days. 
We have shown that the pulse profile has a sinusoidal-like single-peaked shape and 
has practically no dependence on the source luminosity or energy band. 
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INTRODUCTION 



It is conventionally believed that typical X-ray binary pul- 
sars with Be optical companions (hereafter Be/XRP) mani- 
fest themselves through transient activity of two types (see, 
e.g., H (loll])): 

- Type I outbursts are caused by the increase of the 
mass accretion rate onto the neutron star during the perias- 
tron passage. Such flares are usually characterized by their 
periodic appearance (once per binary orbit), short duration 
(small part of the orbital period) and peak X-ray luminosity 
about or less than ~ 10 37 erg s _1 ; 

- Type II outbursts deal with the non-stationary in- 
crease of an amount of matter in the circumstellar disc 
around the Be star and can occur at any orbital phase. The 
duration of such events varies from weeks to months, during 
which the source X-ray luminosity can reach the Eddington 
limi t (~ 10 38 erg s" 1 )- 

iReig fc Rochel l| 19991 ) pointed out the existence of a sub- 
class of binaries with Be companions which harbor a slowly 
rotating neutron star and are characterized by a persistent 
low-luminosity X-ray emission. Those authors proposed also 
the model for such sources in which the neutron star. orbit- 
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ing around a Be star in a relatively wide binary system, 
accretes material only from the low-density outer regions of 
the normal star envelope. Studies of such long-period low- 
luminosity sources are crucial to understand the Be/XRP 
binary systems formation processes and investigations of the 
interaction between the Be star circumstellar disc and the 
neutron star. 

The X-ray source RX J0440. 9+4431 was found dur- 
ing the ROSAT Galactic plane survey with the optical 
companion BSD 24 -491/L S V +44 17 classified as a Be star 
|Motch et al.ll 1997ft . Later lReig fc Rochel (|l999T ) detected X- 
ray pulsations from this source with the period of ~ 202.5 
s and related it to the not numerous family of persistent 
Be/XRP binaries with neutron stars. Based on the optical 
observations t he distance to t he system was estimated as 
3.3 ± 0.5 kpc jReig et alj|2005h . 

The first evide nce of the pu l sar ou tburst activity in X- 
rays was found by iMorii et all (|201Ch in late March 2010 
with the MAXI all-sky monitor. The next outburst was re- 
peated af ter ~ 5 months of a quiescence state in Septem- 
ber 2010 (|Krivonos et alj|2010ah . During the latter flare the 
pulsar was in the INTEGRAL observatory field of view, 
that gave us the opportunity to investigate the properties of 
RX J0440.9+4431 in hard X-rays for the first time. The third 
outburst detected by the Swift observatory in late January 
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Figure 1. One-day averaged light curve of the X-ray pulsar 
RX J0440.9+4431 obtained with the Swift/BAT monitor in the 
15-50 keV energy band (top panel). The thin upper marks corre- 
spond to the RXTE observations, the thick lower marks - to the 
Swift/XRT observations. The light curve obtained by the INTE- 
GRAL/IBIS telescope in the 20-60 keV energy band is shown in 
the inset as a comparison to the Swift/BAT points. 



2011 allowed iTsygankov et all (|201ll ) to propose the exis- 
tence of the ~ 155-days orbital period in the binary system. 
During these outbursts the X-ray luminosity of the source 
reached ~ 9 x 10 36 , ~ 3 x 10 36 and ~ 2 x 10 36 erg s -1 , 
respectively (hereafter "luminosity" means the luminosity 
in the 3-100 keV energy band). These values are about two 
orders of magnitude higher than the quiescent luminosity 
~ 4 x 10 34 erg s _1 , see Fig(TJ . This formally moves formally 
RXJ0440.9+4431 fro m low-lumi nosity persistent Be/XRB 
class to transient one (|Reidl201ll ). 

In this paper, based on the data of Swift, RXTE and 
INTEGRAL observatories in a wide energy band (0.6 — 120 
keV), we report the discovery of a cyclotron absorption line 
in the spectrum of the X-ray pulsar RX J0440. 9+4431 at 
the energy ~ 32 keV, that corresponds to a magnetic field 
strength of 3.2 x 10 12 G (Section 3). A comprehensive analy- 
sis of variations of the source power density spectrum (PDS) 
with the luminosity during type I outbursts in April and 
September 2010 confirmed this estimation (Section 4.1). Re- 
sults of the timing analysis on different time scales (pulse 
and orbital periods) are presented in Sections 4.2 and 4.3. 



INTEGRAL/ISGRI flux measurements in the 20 - 60 keV 
energy band are represented by triangles in the inset. 

Observations in the soft energy band 0.8 — 9 keV were 
performed with the XRT teles cope onboard the Swift Ob- 
servatory (Ge hrels et alj|2004t ) in different intensity states 
(Obs.IDs 00031690, 00418109 and 00418178, thick lower 
strokes in Fig[TJ both in the Photon Counting (PC) and 
Window Timing (WT) modes. As the observations in the 
PC mode performed in a high-luminosity state were affected 
by the pile-up, in the subsequent analysis we used only data 
obtained in the WT mode for uniformity. The typical ex- 
posure during XRT observations was about of 3 — 5 ksec 
each. 

In the standard X-ray energy b and (4 - 20 keVU data 
from the RXTE/PCA spectrometer (|Bradt et al.ll 19931 ) were 
used (Obs.ID 95418-01-XX-XX). The spectral analysis was 
done using the Standard-2 data, temporal analysis - using 
the Good Xenon data with a high time resolution. Note 
that all RXTE observations cover mostly the fading phase 
of both outbursts (thin upper strokes in Fig[T} and have an 
exposure of several kiloseconds. 

This work also employs the INTEGRAL observatory 
(|Winkler et al.l 120031 ) data acquired during the Galactic 
latitude scan campaign at 1=155° (Proposal ID 0720049, 
PI R. Krivonos). The source RX J0440. 9+4431 was in the 
field of view of both X-ray telescopes of the observatory: 
JEM-X (energy band 6-20 keV) and IBIS/ISGRI (energy 
band 18 - 120 keV). The reduction of the ISGRI detector 
da ta was done using recen tly developed methods described 
by iKrivonos et all |2010bl ). The source spectrum from the 
JEM-X telescope was extracted with the standard OSA 
package version 9.(0. Both spectra (from JEM-X and IBIS 
telescopes) were averaged for several days near the maxi- 
mum of the September 2010 outburst with total exposures 
of ~ 47 and ~ 180 ksec, respectively. Such a difference 
in the exposure time is connected with an observational 
pattern and d i fferent fields of view of the instruments (see 
IWinkler et al. | (|2003l ) for details). The method used for IS- 
GRI data spectra extraction has a systematic error in the 
source flux determination of about 3%. We included this 
systematic uncertainty in spectral analysis done with the 
XSPEC package. 

The final spectral and timing analysis for all in- 
struments were done using the standard tools of the 
FTOOLS /LHEASOFT 6.7 package. 



3 SPECTRAL ANALYSIS 
3.1 The Continuum Spectrum 

The spectrum of RX J0440 9+4431 in the 3 - 2 keV en- 
ergy band was obtained bv lReig fc Roche] l| 19991 ) using the 



2 OBSERVATIONS AND DATA ANALYSIS 

To get the broadband view on the properties of the X-ray 
pulsar RX J0440. 9+4431 we used data from three currently 
operating observatories: INTEGRAL, Swift and RXTE. The 
data collected during both outbursts in April and September 
2010 as well as in a quiescent state were analyzed. The source 
light curve obtained with the Swift/BAT telescope in the 
15 — 50 keV energy band is shown in Fig[T]by circles. The 



1 The latest improvement in the response matrix 
(http://www.univcrsc.nasa.gov/xrays/programs/rxte/pca/doc/ 
rmf/pcarmf-11.7/) allows formally to extend the PCA upper 
energy limit up to 50-60 keV. But our analysis showed that there 
are still significant systematic deviations at the energy of ~ 30 
keV. Therefore, in this paper we limited oneself to using just a 
standard energy band. 

2 http://isdc.unigc.ch 
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Table 1. Best-fit results in the energy range 0.6-120 keV 



Parameter 


Model I 


Model II 


N H X 10 22 , cm" 2 


0.47 ± 0.04 


0.46 ±0.04 


kT BB , keV 


1.57 ±0.03 


1.56 ±0.03 


Rbb, km 


0.91 ±0.04 


0.93 ±0.05 


Photon index 


0.74 ±0.05 


0.75 ±0.05 


E cut , keV 


16.9 ± 1.0 


18.4 ±1.1 






0.37 ±0.06 


E C ycl, keV 




31.9 ± 1.3 


<?cycl> keV 




6 (fixed) 


RFjEM-X 


1.08 ±0.02 


1.08 ±0.02 


RFpca 


0.55 ±0.02 


0.55 ±0.02 


RFlBIS 


1.16 ±0.06 


1.26 ±0.06 


X 2 (d.o.f) 


1.315(331) 


1.223(329) 



All uncertainties correspond to one standard deviation 



Figure 2. Energy spectrum of RX J0440.9±4431 obtained during 
the outburst in September 2010 with XRT (green points), JEM- 
X (magenta points), PC A (blue points) and IBIS (red points) 
telescopes (a). The residuals to the fit in units of sigma without 
(b) and with (c) a ~ 32 keV cyclotron absorption line. 

RXTE observatory data acquired in 1998 in the quiescent 
state (the corresponding luminosity in 3-30 keV energy band 
was ~ 3 x 10 34 erg s" 1 ). The source faintness didn't allow 
those authors to distinguish meaningfully between the dif- 
ferent spectral models or determine unambiguously spectral 
parameters. Particularly, the value of the photoelectric ab- 
sorption was varied from its absence to ~ 6 x 10 22 atom cm" 2 
depending on the model, whereas the interstellar absorption 
based on the measurements of the 21 cm emission in the 
sourc e direction i s ~ 0.6 x 10 22 atom c m" 2 (|Kalberla et alj 
120051 ). Moreover, [keig fc Rochel (|l999l ) also reported about 
a marginal detection of a fluorescent iron line with an upper 
limit to its equivalent width of ~ 100 eV. 

On the other hand proper knowledge of the contin- 
uum is crucial not only to understand the physical mech- 
anisms of the X-ray emission generation but also detect 
possible line-like features in the source spectrum. There- 
fore, we first investigated the spectrum of the continuum 
of RXJ0440. 9+4431. 

As mentioned above the XRT telescope observed 
RX J0440. 9+4431 several times both in the quiescent state 
and during the outbursts. In all cases, independently of the 
spectral model used, the source spectrum didn't require the 
inclusion of a strong photoelectric absorption. Its value var- 
ied insignificantly in the range of Ah — (0.2 — 0.5) x 10 22 
atom cm -2 , which is compatible with the value of interstellar 
absorption in the source direction. The quiescent spectrum 
was well fitted by a powerlaw model with an exponential 
cutoff (cutoffpl model in XSPEC), that is typical for X-ray 
pulsars, but the required cutoff energy was unusually low 
(~ 1.8—2.4 keV). The source spectra obtained in a high state 
could not be described adequately by a single-component 
model. The addition of a black-body component with a tem- 
perature of kT B B — 1.4 — 1.7 keV to the powerlaw model 
significantly improved the fit quality. As the detailed study 
of the spectral evolution of RXJ0440. 9+4431 is beyond the 
scope of this paper, we give here only as an example the 



spectral parameters obtained during the high state in April 
2010 (Obs.ID 00031690001, MJD 55299.9): temperature of 
the black-body component kTBB = 1-50 ± 0.05 keV, photon 
index V = 1.13 + 0.27, absorption A H = (0.46 + 0.07) x 10 22 
atom cm - . 

Widening the considered energy band up to 0.6 — 20 
keV by the inclusion of PCA/RXTE data doesn't change 
the described picture or spectral parameters significantly; 
the only thing that was added is the exponential cutoff with 
energies of 5 — 8 keV. It is important to note that changes 
of the source flux during outbursts by a factor of ~ (2 — 3) 
didn't lead to changes of the shape of the source spectrum, 
therefore we jointly fitted XRT and PCA data just adding 
a normalization factor. 

Finally, we could not detect the fluorescent iron line at 
6.4 keV in any luminosity state and obtained only an upper 
limit to its equivalent width ~ 50 eV (90% confidence level). 



3.2 Detection of a Cyclotron Absorption Line 

The cyclotron resonance scattering features (CRSFs) with 
typical energies of ~ 20 — 50 keV are observed at the moment 
in mo r e than a dozen X-ray p ulsars (see, e.g.. lCoburn et all 
l|2002l ); lFiiippova et all |2005l )). These features are thought 
to be due to photon resonant scattering by electrons at Lan- 
dau orbits in the strong magnetic field (of the order of 10 12 
G) of the neutron star and, thus, their registration in the 
source spectrum gives a direct measurement of the magnetic 
field value. 

The pulsar RX J0440. 9+4431 was detected for the first 
time in the hard X-ray band (> 20 keV) by the IBIS tele- 
scope during the September 2010 outburst. Since the source 
flux was relatively low, only the spectrum averaged over the 
whole outburst had enough statistics for the following anal- 
ysis. To reconstruct the broadband (0.6 — 120 keV) spec- 
trum we used data from XRT (Obs.ID 00031690003, MJD 
55443.03), JEM-X and IBIS (MJD 55440-55449) telescopes; 
an averaged flux in the 3-100 keV energy band during these 
sets of observations was ~ 2.3 x 10~ 9 erg cm" 2 s" 1 . To 
get more strict constraints on the spectrum parameters we 
introduced to the fitting procedure also PCA data (Obs.ID 
95418-01-03-00, MJD 55447.96). It is necessary to note, that 
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this observation was performed several days after the out- 
burst maximum, when the source flux was several times less. 
But, as it was mentioned in §3.1, the spectrum shape re- 
mains practically the same with changes of the source flux, 
therefore we can jointly fitted all these data leaving the nor- 
malization factor free. 

The energy spectrum of RX J0440. 9+4431 is shown in 
FigEK (green points - XRT data, magenta points - JEM-X 
data, blue points - PC A data and red points - IBIS data). 
Based on the previous analysis (§3.1) it was approximated by 
the two-component model (black-body radiation and pow- 
er law with an exponential cutoff) with a photoelectric ab- 
sorption (Model I). Best-fit parameters are listed in Ta- 
ble[T] (RFjeai~x, RFpca and RFibis are renormalization 
factors between different instruments, RFxrt = 1.0, was 
fixed). Note, that a difference between the obtained photon 
index and the cutoff energy from the best-fit results in the 
standard energy band is due to the inclusion of higher energy 
data. This model describes the data relatively well every- 
where except around ~ 30 keV, where one can clearly see an 
absorption like deviation from the continuum model (Fig(2j 
middle panel). It is worth noting that the appearance of 
this feature is not a sequence of the model of the continuum 
used. We investigated several other combinations of spec- 
tral components to describe the continuum (NPEX + black- 
body, highecut+ black-body) and found that: 1) all of them 
approximate it as well as cutoffpl+ black-body, but have 
more parameters; 2) residuals for all of them demonstrate 
a prominent absorption- like feature near 30 keV. We also 
checked that this feature is not an artifact connected with 
the response matrix of the ISGRI detector. For this purpose 
spectra of a number of sources with a similar intensity and 
without spectral features in the range of interest were re- 
constructed from nearby observations of the Galactic Bulge. 
No unusual absorption features were found in the spectra 
of these sources. Thus, we can conclude that the absorption 
feature in the spectrum of RX J0440. 9+4431 is the attribute 
of the source arising due to a cyclotron resonant scattering. 

To describe the detected feature three different spectral 
models were used: 

1. multiplicative XSPEC model cyclabs in the form of: 



exp 



ol) 2 +<r 



5 

cycl 



(1) 



where E cyc i, cr cyc i, and r cyc i are the line c e ntral energy, 
width, and depth, respectively (|Mihara et all (|l990T l); 
2. multiplicative XSPEC model gabs in the form of: 



(E-E aycl y 

cycl 



(2) 



(see eq. 6 and 7 in lCoburn et all (20Q2)); 

3. simple additive Gaussian with a negative normaliza- 
tion. 

Best-fit parameters for the case of the cyclabs model 
(Model II) are summarized in Table Q] It is seen that an 
addition of the cyclotron absorption line leads to a sub- 
stantial reduction by \ 2 — 30 for 2 d.o.f. Nevertheless, to 
estimate the significance of the line detection we performed 
10 Monte-Carlo simulations of the spectra based on the 
best-fit model without a line. Then we fitted the obtained 
spectra adding the cyclabs component and analyzed the dis- 
tribution of the line depth (a similar approach was used by 



iRea et all (|2005h ). Parameters of the model for the simu- 
lated spectra were the same as for the real spectrum: the 
line width was fixed at 6 keV, but the line depth and energy 
were free parameters. None of the absorption lines detected 
in the simulated spectra was deep as in the observed one. 
Thus, based on these simulations we can conclude that the 
significance of the cyclotron absorption line in the spectrum 
of RX J0440. 9+4431 is higher than ~ 4a. 

The best-fit positions of the line energy in other cases 
depend on the model used, but nevertheless agree with the 
cyclabs results within error bars: E cyc i = 34.2 ± 1.5 keV for 
the gabs model and E cyc i = 30.6 ± 1.6 keV for the Gaussian 
one. It is important to note that in all three models we 
could not correctly constrain the line width, probably due 
to a source faintness and small depth of the line, therefore it 
was fixed at the value of 6 keV, that is typica l for cyclotron 
lines i n the spectra of X-ray p ulsars (see, e.g.. lCoburn et all 
(|2002l ), lFilippova et all (|2005h ). 

From the measured value of the cyclotron line energy we 
can estimate the magnetic field on the neutron star surface 
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3.2 x 10 12 G 



where Rns = 15 km ()Suleimanov et aI1l201ll ) and Mj\ 



I AM, 



© 



are the neutron star radius and mass estimates, 



respectively. 

Finally, it is necessary to note that despite the signifi- 
cant improvement of the fit quality an excess of emission at 
energies higher than ~100 keV still remains. Formally it can 
be fitted by an additional powerlaw component, but due to 
the lack of statistics at such high energies it is impossible to 
obtain significant restrictions to its parameters and to make 
conclusions about its nature. 



4 TIMING ANALYSIS 
4.1 Power Spectra 

Accreting X-ray pulsars are characterized by a specific shape 
of the noise power spectrum. Due to the presence of the 
strong magnetic field of the neutron star, the disk-like ac- 
cretion flow is disrupted at the magnetospheric boundary. 
This divides the flow into two distinct parts - the accretion 
disk and th e magnetospheric flow, which have different noise 
properties (|Revnivtsev et aL 20091) . 

It was shown bv lRevnivtsev et al.l ([20091 ) that the break 
frequency in the noise power spectra of accreting X-ray pul- 
sars reflects the timescale of noise generation at the inner 
boundary of the accretion disk/flow, and its value depends 
on the mass accretion rate in the binary system. An increase 
of the mass accretion rate reduces the size of the magneto- 
sphere (and hence the inner radius of the disk), so that the 
characteristic frequency at the inner edge of the disk/flow 
is increased. This property of the power spectrum can be 
used for estimations of the magnetic moment of the accret- 
ing compact star. 

Due to the source faintness and short duration of the 
observations it was difficult to build the qualitative power 
spectrum for each RXTE observation. But we succeeded 
by averaging them into four points with mean luminosity 
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Figure 3. Top: Power density spectra of RXJ0440. 9+4431 in dif- 
ferent luminosity states. The spectrum in the bright state is mul- 
tiplied by a factor of 50 for clarity. The frequency shown along the 
X-axis is expressed in units of the compact object spin frequency. 
The power density spectra are multiplied by the frequency. The 
solid lines show the best-fit models (see text for details). Bottom: 
dependence of the break frequency in the noise power spectrum 
of RX J0440. 9+4431 on the bolometric luminosity of the source. 
The dashed line shows the prediction of the model described in 
the text. The dotted line shows the neutron star spin frequency. 



~ 4 x 10 34 erg s" 1 (ID 30111-01-01-00), ~ 1.3 x 10 35 erg 
s" 1 (IDs 95418-01-04-00, 95418-01-04-01), ~ 2.4 x 10 35 erg 
s" 1 (IDs 95418-01-02-01, 95418-01-03-00, 95418-01-03-01) 
and ~ 7.1 X 10 36 erg s" 1 (IDs 95418-01-01-00, 95418-01- 

02- 00) (here we took into account the fact that the spec- 
trum shape is nearly constant with the luminosity and re- 
calculated the PCA flux in the 3-20 keV energy band to the 

3- 100 keV one). In Fig[3] (top panel) two power spectra of 
RXJ0440. 9+4431 are shown for different luminosities. The 
blue histogram corresponds to the PCA observation 30111- 
01-01-00 (MJD 50843.27) in a low luminosity state with 
Lx — 4 x 10 34 erg s - , the red histogram represents the 
bright state with an averaged luminosity of Lx ~ 7.1 x 10 36 
erg s _1 (multiplied by a factor of 50 for clarity). The fre- 
quency shown along the X-axis is expressed in units of the 
compact object spin frequency: the pulse period was 202.5 
l|Reig fc Rochdl 19991 ) and 205.0 s (this work) in the low and 
high luminosity s t ates, r espectively; see Sect. 4.2 for details. 
iRevnivtsev et al.l (|2009D showed that the noise power spec- 
tra of all examined sources have a similar power-law slope 
(P ~ f~ 2 ) above the break frequency, irrespective of the 



power spectra form at lower frequencies. To fit the PDSs of 
X-ray pulsars these authors proposed a simple model in the 
form of P oc / X (1 + (///b) 4 )^ ' 5 . 

Solid lines in Fig[3] represent results of the approxi- 
mation of the power spectra of RX J0440. 9+4431 with the 
described model: in the low-intensity state one-component 
model was used; in the high-intensity state the second com- 
ponent in the same form was added to the model. Thereby 
an additional noise component was taken into account. This 
component is presumably generated in the ring of the accre- 
tion disk between the radii corresponding to the size of the 
magnetosphere at the high accret ion rate (small radius) a nd 
low accretion rate (large radius) l|Revnivtsev et al]|2009l ). 

Taking the expression for the frequency of the Keplerian 
rotation as 2-kvk = (GM) 1 ^ 2 R^^ 2 and for magnetospheric 
radiu s as R M = (2GM)~ 1/7 M~ 2/r (|Bildsten et al.1 

Il997t ) one can show that the break frequency then depends 
on the neutron star's magnetic moment and mass accretion 
rate as: 

/fc = ^(GAf) 5/ V 6/7 M 3 / 7 , (3) 

where ft is the break frequency, /j, is the dipole magnetic 
moment of the neutron star and M is the mass accretion 
rate. 

Substituting in Eq[3] R NS = 15 km, M NS = 1.4M Q , 
3.3 kpc for the distance to the source and L x = 0.1c 2 M we 
can approximate the measured values of the break frequency 
(dashed line in the bottom panel of Fig[3]) and determine 
the strength of the magnetic field of the neutron star as 
B = 3.3 x 10 12 G, which is in very good agreement with 
the one derived from the cyclotron absorption line energy 
detected in the source spectrum. 

4.2 Pulse Profile 

It is well known that the pulse profiles of X-ray pulsars are 
stro ngly dependent on photon ener gy and luminosity (see, 
e.g.. lLutovinov fc Tsygankovl (|2009l ) and references therein). 
This is especially prominent in a soft energy band strongly 
affected by the absorption in the vicinity of the neutron star. 
Therefore the evolution of the pulse profile can provide an 
important information for understanding the properties and 
spatial distribution of matter around the neutron star. To 
reconstruct pulse profiles we determined the pulse period 
during each observation. We didn't detect any significant 
variations of the pulse period during the outbursts therefore 
in the following analysis we used the value averaged over the 
whole September 2010 outburst 205.0 + 0.1 s. This period is 
significantly longer t han the one measured in 1998 ~ 202.5 s 
(|Reig fe Roche|[l999l ), that indicates the deceleration of the 
neutron star rotation over the last decade with an average 
rate of P/P ~ 10 -3 yr -1 . 

In Fig|4]the pulse profile of RX J0440. 9+4431 in dif- 
ferent energy bands is shown to illustrate the evolution 
of it's shape with the energy in the low luminosity state 
(Lx ~ 1.5 x 10 36 erg s _1 ). In general the pulse profile 
of RX J0440. 9+4431 has a sinusoidally-like single-peaked 
shape in a wide energy band. But on the light curves with 
high statistics (obtained with the PCA spectrometer in 3 — 8 
and 8 — 14 keV energy bands) some substructures become 
apparent. Particularly, one can see at least two prominent 
features: subpulses before and after the main peak. 
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Figure 4. Pulse profiles of RX J0440. 9+4431 in different energy 
bands (normalized by a mean flux). Data from Swift/XRT (obser- 
vation 00031690003, MJD 55443.03), RXTE/PCA (observation 
95418-01-04-00, MJD 55449.04) and INTEGRAL /ISGRI (revolu- 
tions 963, 964, 965; MJD 55441.1-55448.4) were used (an average 
luminosity is about 1.5 X 10 36 erg s _1 ). The profile is shown twice 
for a clarity. 



In the high luminosity state (the mean luminosity is 
Lx — 8 x 10 36 erg s _1 ) the pulse profile has approximately 
the same shape, but shows another feature - a dip-like struc- 
ture at a phase of ~ 0.25 after the main peak. Th is feature is 
clearly seen at t he energies below 8 keV (Fig|5J.|u sui et al.l 
(|in preparation! ) interpreted this feature as a partial eclipse 
of the emission region by an accretion column of the neu- 
tron star. It is interesting to note that in the low luminosity 
state (Fig(4} such an absorption feature wasn't detected in 
the pulse profile. Probably its appearance can be connected 
with changes of the geometry of the accretion column or 
with changes of the neutron star orientation relative to the 
observer. 
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Figure 5. The same as in FigE] but for the high-luminosity state 
with L x =i 8 X 10 36 erg s _1 . Data from Swift/XRT (observa- 
tion 00418178000, MJD 55289.85) and RXTE/PCA (observation 
95418-01-01-00, MJD 55292.33) were used. The profile is shown 
twice for a clarity. 

The pulsed fractioijf] of the source emission in the stan- 
dard X-ray band (4-20 keV) is relatively high: (40.9±0.1)% 
at high luminosities (~ 8 x 10 36 erg s~\ MJD 55292.33) and 
increases with the fading of the pulsar intensity. In partic- 
ular, in observations 95418-01-04-00 (MJD 55449.04) and 
95418-01-04-01 (MJD 55450.23), when the source luminos- 
ity dropped to ~ 1.3 x 10 36 erg s _1 and ~ 9 x 10 35 erg s~\ 
the pulsed fraction came up to the values of (50.1±0.1)% 
and (61.5±1.0)%, respectively. Such a behavior is typical 
for X-ray pulsars (|Lutovinov fc Tsygankoy||2009l ). However 
for RX J0440. 9+4431 the pulsed fraction shows a tendency 
to decrease towards higher energies (FigJBJ. 



4.3 Orbital Period 

In 2010-2011 three outbursts from RX J0440. 9+4431 were 
sequentially detected by the Swift observatory (see FigJlJ. 
The peak luminosities of all three were relatively low: ~ 



3 It was determined as PF = (7 max — I-mia) / (Jmxx + I-mia), where 
7 max and / m i n are maximum and minimum intensities in the pulse 
profile of an X-ray pulsar, respectively. 
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Figure 6. Pulsed fraction dependence on the energy for two in- 
tensity states: L x — 1.3 X 10 36 erg s _1 (open circles, RXTE 
observation 95418-01-04-00, MJD 55449.04) and L x ~ 8 x 10 36 
erg s" 1 (filled circles, an RXTE observation 95418-01-01-00, MJD 
55292.33). 



9 x 10 36 , ~ 3 x 10 36 and ~ 2 x 10 36 erg s" 1 (assuming 
distance to the source of 3.3 kpc), that allows us to con- 
sider them as Type I outbursts. It was mentioned above 
that such events are separated by one orbital period of the 
binary system. Taking into account that all three outbursts 
were roughly equally spaced in time we estimated the orbital 
period of RXJ0440. 9+4431 as ~ 155 days (a mean time 
between starts of the outbursts). This value is well agreed 
with the period of about 150 d ays derived fro m the Corbet 
diagram of P sp i„ versus P or ut ((Co rbet 1986;). An accurate 
determination of the orbital parameters of the binary sys- 
tem should be based on measurements of the Doppler delay 
of the pulse period due to the orbital motion of the neutron 
star in the binary system and requires much longer observa- 
tions in a combination with the high sensitivity and timing 
accuracy. 



5 CONCLUSION 

In this work we have presented the temporal and 
spectral analysis of the poorly studied X-ray pulsar 
RXJ0440. 9+4431 with the Be star companion. Using the 
data of Swift, RXTE and INTEGRAL Observatories the 
broadband spectrum of the source was reconstructed for 
the first time. It was shown that the spectrum can be well 
approximated by the two-component model, including the 
black-body emission with the temperature of &Tbb — 1.5 
keV and the high-energy component in the form of a power- 
law with an exponential cutoff. Such a spectral shape is sim- 
ilar to that observed for other low-luminosity X-ray pulsars 
in Be binary systems, like X Persei or RX J0146. 9+6121, 
and reminiscent of the thermal emission from the neu- 
tron star polar caps and the co mptonized non t herma l 
emission at high energies (see, e.g ., Idi Salvo et al.l (|l998l ); 
lLa Palombara fc Mereghettil (|2006l ) and references therein). 
The corresponding black-body emitting radius Rbb (see Ta- 



ble!]]) is about 1.1 km which agrees with the expected radius 
of an acc retion column of t h e neu tron star R co i ~ 0.1i?jvs 
(see, e.g. jBasko fc Sunvaevl (|l976t )). 

We have discovered a prominent absorption feature at 
the energy ~ 32 keV in the spectrum of the source. Inter- 
preting it as a cyclotron resonance scattering feature the 
magnetic field strength of the neutron star was estimated 
as B ~ 3.2 x 10 12 G. The obtained value is in very close 
agreement with one derived from the properties of the noise 
power spectra. 

The pulse profile of RX J0440. 9+4431 has a 
sinusoidally-like single-peaked shape, which is stable 
in relation to the changes of the source luminosity and 
energy band. Nevertheless in the high intensity state we 
found a dip-like structure at the pulse phase of ~ 0.25 
after the main peak that can be interpreted as a partial 
eclipse of the emission region by the accretion column of 
the neutron star. The pulsed fraction of the source emission 
is growing with the decreas e of the source luminosity, tha t 
is typical for X-ray pulsars (|Lutovinov fc Tsygankovll2009l ). 
On the other hand, in contrast to other X-ray pulsars, for 
RX J0440. 9+4431 the pulsed fraction shows a tendency to 
decrease towards higher energies. 

Based on the recurrence time between Type I outbursts 
the orbital period of ~ 155 days was supposed for the bi- 
nary system. The gradual decrease of the total energy release 
during these transient events can be probably connected 
with the existence of the circumstellar disk around the B e 
star and its va r iabilit y in time (see, e.g.. IClark et all (|2003l ); 
iRivinius et all (|200ll )). 
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